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Abstract 

The  seasonal  variability  of  the  Japan/East  Sea  (JES)  circulation  was  studied  numerically 
using  the  Princeton  Ocean  Model  (POM)  with  horizontal  resolution  varying  from  1  I  54  to  18.53 
km  and  15  sigma  levels  conforming  to  a  relatively  realistic  bottom  topography  The  model  was 
integrated  using  climatological  monthly  mean  wind  stresses,  heat  and  salt  fluxes  as  surface  forcing 
and  observational  oceanic  inflow/outflow  at  open  boundaries,  The  seasonally  averaged  effects  of 
isolated  forcing  terms  are  presented  and  analyzed  from  the  following  experiments:  1)  non- linear 
effects  removed.  2)  no  lateral  transport  at  open  boundary,  and  3)  wind  effects  removed  Major 
currents  are  simulated  reasonably  well  compared  to  observations.  The  nonlinear  advection  does 
not  affect  the  general  circulation  pattern  evidently,  but  does  affect  the  formation  of  the  mesoscale 
eddies,  especially  the  Ulleung/Tsushima  Dasin  (UTB)  eddy  (all  seasons)  and  the  Japan  Basin 
(JB)  cyclonic  gyre  (spring)  The  lateral  boundary  forcing  enhances  (weakens)  the  JES  volume 
transport  in  the  summer  (winter),  The  wind  forcing  is  the  most  important  factor  (SOSc)  for 
generating  the  JB  cyclonic  gyre  Besides,  it  drives  tile  Liman  Current  and  damps  the  EasI 
Korean  Warm  Current  in  the  winter,  and  generates  the  UTB  eddy,  and  eddies  along  the  Japan 
Nearshore  Branch  (JNB)  in  all  seasons,  However,  it  has  almost  no  effect  on  the  JNB  currents 
for  all  seasons 

1  Introduction 

The  Japan  Sea,  known  as  the  East  Sea  in  Korea,  has  a  steep  bottom  topography  (Fig. 
1 ;  that  makes  it  a  unique  semi-enclosed  ocean  basin  overlaid  by  a  pronounced  monsoon 
surface  wind.  The  Japan/East  Sea,  hereafter  referred  to  as  JES,  covers  an  area  of  10® 
km2.  H  has  a  maximum  depth  in  excess  of  3,700  m,  and  is  isolated  from  open  oceans 
except  for  small  (narrow  and  shallow)  straits,  which  connects  the  JES  witli  the  North 
Pacific  through  the  Tsushima/Korean  and  Tsugaru  Straits  and  with  the  Okhotsk  Sea 
through  the  Soya  and  Tatar  Straits.  The  JES  contains  three  major  basins  called  the 
Japan  Basin  (JB),  Ulleng/Tsushima  Basin  (UTB),  and  Yamato  Basins  (YB),  and  a  high 
central  seamount  called  the  Yamato  Rise  (YR).  The  JES  has  great  scientific  interest 
as  a  miniature  prototype  ocean.  Its  basin-wide  circulation  pattern,  boundary  currents, 
Subpolar  Front  (SPF),  mesoscale  eddy  activities  and  deep  water  formation  are  similar  to 
those  in  a  large  ocean. 

Numerical  studies  on  the  JES  circulat ions  started  in  early  1980.  Various  types  of  mod¬ 
els  were  used  such  as  the  multi-layer  model  (Sekine,  1986,  1991;  Kawabe,  1982b;  Yoon, 
1982a, b;  Seung  and  Nam,  1992a. h:  Seung  and  Kim,  1995),  the  Modular  Ocean  Model 
(MOM)  (Kim  and  Yoon,  1994;  Holloway  et  al.,  1995;  Kim  and  Yoon,  1996:  Yuehikawa  ct 
al.,  1999),  rigid-lid  z-level  model  (Yoshikawa  et  al  ,  1999),  the  Miami  Isopycnul  Coordinate 
Model  (MICOM)  (Seung  and  Kim,  1993)  and  the  Princeton  Ocean  Model  (POM)  (Mooer 
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and  Kang,  1995;  Chu  <■(  al .  1999c,d).  Mont  of  the  numerical  efforts  nre  concentrated  on 
simulating  basin-wide  circulation,  the  Tsushima  Warm  Current  (TWC)  bifurcation,  and 
formation  of  the  intermediate  waters.  However,  there  is  lack  of  studies  on  dynamical 
mechanisms  for  the  formation  of  seasonal  variability  of  the  JES  circulation 

In  thfc>  study,  we  use  the  POM  to  inviistigate  driving  mechanisms  for  the  formation 
of  the  seasonal  variabilities  of  the  IES  circulation  including  the  subpolar  front  (SPF) 
meandering  and  eddies,  the  TWC  bifurcation  and  its  effect  on  the  formation  of  niwsocif 
eddies  in  the  1TB  and  the  YB,  and  the  Liman  Cold  Current  (LCC)  and  its  penet  rat  ion 
into  the  southwestern  waters  along  the  Korean  coast  The  control  run,  forced  by  the  cli¬ 
matological  monthly  wind  stress,  heat  and  fresh  water  fluxes,  is  designed  to  best  simulate 
reality  against  which  each  experiment  is  compared.  In  the  experiments,  various  external 
and  internal  fac  tors  nre  modified  and  the  resulting  circulation  patterns  and  magnitudes 
compared  to  the  control  run  results.  Spec  ifically,  we  estimate-  the  contribution,  in  terms 
of  volume  transport  and  circulation  patterns,  of  non-linear  advection,  wind  forcing  and 
Intend  boundary  transport  to  the  ocean  features  identified  in  the  control  results  From 
this,  wc  can  estimate  the  relative  importance  of  these  factors  to  the-  seasonal  variability 
of  the  JES  circulation 

2  JES  Current  Systems 

Must  of  the  nearly  homogeneous  water  in  the  deep  part  of  the  basin  is  called  the  Japan 
Sea  Proper  Water  ( Moriya.su,  1972)  and  Is  of  low  temperature  and  low  salinity.  Aliove 
the  Proper  Water,  the-  TWC,  dominating  the  surface  layer,  flows  from  the  East  China  Sea 
through  the  Tsushima/Korean  Strait  and  carries.  warm  water  from  the  south  The  LCC 
flows  in  the  JES  from  the  Okhotsk  Sea  through  the?  Tatar  Strait,  carries  cold  fresh  water 
from  Ihc  north,  and  becomes  the  North  Korean  Cold  Current  (NKCC)  after  reaching  the 
North  Korean  coast  (Yoon.  1982)  Both  currents  turn  eastward  to  flow  roughly  along 
tt)  N  latitude,  forming  the  SPF  between  the  Tsushma  Warm  Water  and  the  cold  and 
fresh  water  from  the  north  (Fig.  2). 

Uda  ( 19.11)  was  the  first  one  to  sketch  the  JES  general  circulation  from  limited  oliser- 
vational  data,  The  TWC  separates  north  of  -15  N  into  two  branch**  into  a  western  and 
an  eastern  channel  (Kawabe.  I982a,b;  Hose  el  al..  1999)  and  flows  through  the  western 
channel,  called  the  East  Korea  Warm  Current  (KKWC).  and  closely  follows  the  Korean 
coast  until  it  separatist  near  .'17  N  into  two  branch**.  The  eastern  branch  follows  the  SP1- 
to  the  western  coast  of  Sapporo  Island,  and  the  western  branch  moves  northward  and 
forms  a  cyclonic  eddy  at  the  Eastern  Korean  Bay  (EKB).  It  flows  through  the  eastern 
channel  which  closely  follows  the  Japanese  Coast,  called  the  Nearshore  Branch  (NB)  by 
Yoon  (1982a).  More  accurately,  we  call  it  the  Japan  Nearshore  Branch  (JNB),  The  JNB  is 
usually  weaker  than  the  EKWC.  Tltc  strength  of  the  Tsushima  Current  at  both  channels 
reduces  with  depth. 

The  NKCC  meets  the  EKWC  at  about  .'18’  N  with  some  seasonal  meridional  migration 
After  separation  from  the  coast.  tl»e  NKCC  and  the  EKWC  converge  and  form  a  strong 
front  that  stretches  in  a  west-east  direction  across  the  basin.  The  NKCC  makes  a  cyclonic 
recirculation  gyre  in  the  north  but  most  of  the  EKWC'  flows  out  through  the  outlets. 
I  he  formation  of  NKCC  and  separation  of  EKWC  are  due  to  local  forcing  by  wind  and 
buoyancy  flux  (Seung,  1992)  Large  meanders  develop  along  the  front  and  arc  rsmh  iated 
with  warm  and  cold  eddies  Readers  may  find  qualitative  depiction  from  a  textbook 
written  by  Tomczak  and  Godfrey  (1994). 

Chu  et  al  (1999b)  identified  tiuijor  feature-  of  the  I hrec-dimcnsinnal  circulation  and 
t lie  volume  transport  from  the  Navy’s  0..V  xO..V  global  monthly  climatological  tempera¬ 
ture  and  salinity  data  set  using  the  P-vector  method  (Chu,  1995;  Chu  et  al.,  1998b'  The 
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transport  pattern  is  largely  determined  by  the  upper  layer  circulation  and  characterized 
by  a  large-scale  cyclonic  recirculation  gyre,  in  which  the  EKWC  and  the  JNB  take  part,  as 
the  inflow-outflow  system,  and  the  NKC’C  in  the  North  At  a  few  hundred  kilometers  off 
the  separation  area,  the  EKWC  makes  an  anticyclonic  gyre.  The  gyre  baomes  stronger 
as  the  EKWC  develops  On  the  other  hand,  the  northern  cyclonic  gyre  is  very  deep  and 
is  most  significantly  in  the  winter  strengthened  by  the  wind  and  buoyancy  flux.  The 
gyre,  or  the  southward  coastal  current  related  to  it,  is  deep  enough  to  intrude  southward 
beneath  the  EKWC  most  of  the  time,  Seung  also  confirmed  the  summertime  presence  of 
counter-current  beneath  the  JNB.  North  of  the  SPF  there  exists  n  cyclonic  gyre  in  the 
JB  usually  called  the  JB  gyre. 

3  Seasonal  Variation  of  Atmospheric  Forcing 

The  Asian  monsoon  strongly  affects  the  thermal  struc  ture  of  the  ,1ES  (Chu  ct  al..  199Ka). 
During  the  winter  monsoon  season,  a  very  cold  northwest  wind  blows  over  the  JES  as  a 
result  of  the  Siberian  High  Pressure  System.  By  late  April,  numerous  frontally-generated 
events  occur  making  late  April  and  May  highly  variable  in  terms  of  wind  specs  Is  and 
number  of  clouds.  During  this  period  storms  originating  in  Mongolin  may  cause  strong, 
warm  westerlies.  By  late-  May  and  early  June,  the  summer  surface  atmospheric  low 
pressure  system  begins  to  form  over  Asia.  Initially  this  low  pressure  system  is  centered 
north  of  the  Yellow  Sea  producing  westerly  winds.  In  late  June,  this  low  begins  to  migrate 
to  the  west  setting  up  the  southwest  monsoon  that  dominate.’  the  summer  months,  The 
winds  remain  variable  through  June  until  Manchurian  low  pressure  system  strengthens. 
Despite  the-  very  active  weather  systems,  the  mean  surface  wind  sjx'cd  over  the  JES  in 
summer  is  between  3  and  I  m/s,  which  is  weaker  than  in  winter.  By  July,  however,  high 
pressure  (the  Bonin  High)  to  the  south  and  the  low  pressure  over  Manchuria  produce 
southerly  winds  carrying  warm,  moist  air  over  the  East  China  Sca/Yellow  Sea  In  the 
summer,  warm  air  and  strong  downward  net  radiation  stabilize  the  upper  layer  of  the 
water  and  causes  the  surface  mixed  layer  to  shoal,  October  is  the  beginning  of  t  he 
transition  back  to  winter  conditions.  The  southerly  winds  weaken  and  let  the  sea  surface 
slope  reestablish  the  winter  pattern. 

The  datasets  used  were  the  objectively  analyzed  fields  of  surface  marine  climatology 
and  anomalies  of  fluxes  of  heat,  momentum,  and  fresh  water.  The  fields  are  derived  from 
individual  observations  in  the  Comprehensive  Ocean-At rnosphere  Data  Set  (COADS) 
from  1 9-IT>  to  December  1989  and  are  analyzed  on  a  I  by  l  grid  (da  Silva  ct  al.,  1991), 

4  The  Numerical  Ocean  Model 

4.1  Model  Description 

Coastal  oceans  and  semi-enclosed  setts  are  marked  by  extremely  high  spatial  and  tempo¬ 
ral  variability  that  challenge  the  existing  predictive  capabilities  of  numerical  simulations 
The  POM  is  a  time  dependent,  primitive  equation  circulation  model  on  a  three  dimen¬ 
sional  grid  that  include  realistic  topography  and  a  free  surface  (Blumberg  and  Mellor, 
1987).  River  outflow  is  also  not  included.  However,  the  seasonal  variation  in  sea  sur¬ 
face  height,  temperature,  salinity,  circulation  and  transport  arc  well  represented  by  the 
model  From  a  series  of  numerical  experiments,  the  qualitative  and  quantitative  effects 
of  non-linearity,  wind  forcing  and  lateral  boundary  transport  on  the  JES  are  analyzed, 
yielding  considerable  insight  into  the  external  factors  affecting  the  region  oceanography. 
The  model  results  were  sampled  every  thirty  days 
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The  model  contains  94  x  100  x  15  horizontally  fixed  grid  points.  The  horizontal  spacing 
js  10’  latitude  and  longitude  (approximately  11.54  to  15.18  km  in  the  zonal  direction 
and  18.53  km  in  the  latitudinal  direction)  and  15  vertical  sigma  coordinate  levels  are 
used.  The  model  domain  is  from  35.0°  N  to  51.5°  N,  and  from  127.0°  E  to  142.5° 
E.  The  bottom  topography  is  the  smoothed  data  from  the  Naval  Oceanographic  Office 
Digital  Bathymetry  Data  Base  5  minute  by  5  minute  resolution  (DBDB5).  The  horizontal 
diffusivities  are  modeled  using  the  Smagorinsky  (1963)  form  with  the  coefficient  chosen 
ro  be  0.2  for  this  application.  The  bottom  stress  Th  is  assumed  to  follow  a  quadratic 
law 

ffc  =  Po^’dIVj.IVj,  (1) 

where  pQ  (  =  1025  kg/m3)  is  the  characteristic  density  of  the  sea  water,  V(,  is  the  horizontal 
component  of  the  bottom  velocity,  and  Co  is  the  drag  coefficient,  which  was  determined 
by  the  similarity  theory'- 

4.2  Surface  Forcing  Functions 

The  atmospheric  forcing  for  the  JES  application  of  the  POM  includes  mechanical  and 
thermohaline  forcing.  The  wind  forcing  is  depicted  by 


dudv 

—  =  (tOT-fOw) 


(2) 


where  (u.  v)  and  (r(lr.  my)  are  the  two  components  of  the  water  velocity  and  wind  stress 
vectors,  respectively.  The  wind  stress  at  each  time  step  is  interpolated  from  monthly 
mean  climate  wind  stress  from  GOADS  (1945-1989).  We  interpolated  the  GOADS  wind 
stress  of  the  resolution  of  l°x  1°  onto  the  model  grid  of  the  resolution  of  10'. 


Surface  thermal  forcing  is  depicted  by 


~  +  alC(fioBS  ~  G) 

QQ 

A"?—  =  -OqSfP  -  E)  +  OtiC(SoBS  S) 


(4) 


where  Mobs  and  Sobs  are  the  observed  potential  temperature  and  salinity,  cp  is  the 
specific  heat,  Qh  is  the  net  surface  heat  flux  (downward  positive),  P  is  the  precipitation 
rate,  and  E  is  evaporation  rate.  The  relaxation  coefficient  C  is  the  reciprocal  of  the 
restoring  time  period  for  a  unit  volume  of  water.  The  parameters  (ah  o2)  are  (G.l)-type 
switches:  0|  —  1.  o2  —  0,  would  specify  only  flux  forcing  is  applied;  o,  -  0,  q2  —  1, 
would  specify  that  only  restoring  type  forcing  is  applied. 

Ghu  (1989)  pointed  out  the  importance  of  using  compatible  surface  wind  anil  thermal 
boundary  conditions.  Chu  et  al.  ( 1 998c)  further  found  that  the  restoring  surface  boundary 
condition  does  not  exist  anywhere  in  the  world  ocean.  Therefore,  in  this  study,  the 
surface  thermohaline  forcing  Is  determined  solely  by  the  fltoc  forcing,  that  is,  at  —  1, 
—  ()  in  (3)-(4).  The  mixing  coefficients  ,  K n ,  and  AT.s  were  computed  using  a  level 
two-turbulence  closure  scheme  (Mellor  and  Yamada,  1982). 

4.3  Lateral  Boundary  Forcing 

Closed  lateral  boundaries,  i.e.,  the  modeled  ocean  bordered  by  land,  were  defined  using  a 
free  slip  condition  for  velocity  and  a  zero  gradient  condition  for  temperature  and  salinity. 
No  advective  or  diffusive  heat,  salt  or  velocity  fluxes  occur  through  these  boundaries. 

At  open  boundaries,  the  numerical  gTid  ends  but  the  fluid  motion  is  unrestricted. 
Uncertainty  at  open  boundaries  makes  marginal  sea  modeling  difficult.  Three  approaches. 


98 


ESTUARINE  AND  COASTAL  MODELING 


Mont  h 

Jan 

Feb 

Mar 

■SB 

May 

Jun 
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Id 

Soya 

-0.12 

-0  20 

-0.15 

-0.20 

-0  35 

-0.55 

EjTjl 

•0  1 

Tatar 

o  2 

0.2 

0.2 

0.2 

0.2 

0.2 

Tsugaru 

-0.68 

-0.38 

-0.60 

-0.45 

-0.60 

-1.05 

-1.65 

-1.73 

-1.2 

Tsushima 

0.7 

0.3 

0.6 

0.4 

0.6 

1.2 

2.0 

2.1 

Wtm 

1.1 

Table  1:  Monthly  variation  of  volume  transport  (Sv)  at  the  lateral  open  boundaries  The 
posltlvc/negative  values  mean  inflow /outflow. 


local-type,  inverse-type,  and  nested  liasin/cnastal  modeling,  are  available  for  determining 
the  open  boundary  condition  (Chit  rt  ah.  1997)  Here,  we  take  the  local-type  approach, 
i.e  ,  to  use  the  radiative  boundary  condition  with  specified  volume  transport  When  the 
water  flows  into  the  model  domain,  temperature  and  salinity  at  the  open  tmundarv  are 
likewise  prescribed  from  the  Navy’s  climatological  data.  When  water  Hows  out  of  the 
domain,  the  radiation  condition  was  applied. 

-£(0.S)  +  Cn£(fl,S)  =  0  (5) 

where  the  subscript  is  the  direction  normal  to  the  boundary. 

Warm  water  enters  the  JES  t  hrough  the  Tsushima/Korean  Strait  wit  li  ( he  TWO  from 
the  East  China  Sea.  and  cold  water  enters  the  JES  through  the  Tatar  Strait  with  the 
Liman  Current  from  the  Sea  of  Okhotsk  (Udn,  1934).  The  water  exits  the  IES  through 
t lie  Tsugaru  and  Soya  straits.  We  use  a  more  recent  estimation  of  the  monthly  mean 
volume  transport,  reported  by  Chu  et  a).  (1999b),  through  the  Tsushima/ Korean  Strait 
with  the  annual  average  of  1.3  Sv,  a  muximum  of  2.2  Sv  in  October,  and  a  minimum  of 
0,3  Sv  in  February.  The  total  inflow  transport  through  the  Tatar  and  Tsushima/Korean 
straits  should  be  the  same  as  the  total  outflow  transport  through  the  Tsugaru  and  Soya 
straits.  We  assume  that  7591  of  the  total  inflow  transport  should  flow  out  of  the  JES 
through  the  Tsugaru  Strait,  and  2591  through  the  Sova  Strait.  The  monthly  volume 
transports  at  open  boundaries  arc  listed  in  Table  1. 

4  4  Initial  Conditions  and  Initialization 

The  model  was  integrated  with  all  three  components  of  velocity  (a,  c.  tc)  initially  set  to 
zero,  and  with  temperature  and  salinity  specified  by  interpolating  climatology  data  to 
each  model  grid  point.  The  model  year  consists  of  360  days  (30  days  per  month),  day 
361  corresponds  to  1  January,  It  was  found  that  90  days  were  sufficient  for  the  model 
kinetic  energy  to  reach  quasi-steady  state  under  the  imposed  conditions  In  order  to  first 
capture  the  winter  monsoon  to  summer  monsoon  transition,  the  model  was  started  from 
30  October,  and  run  to  30  January  next  year  for  spin-up.  After  that  day  (30  January), 
the  model  was  run  another  three  years  for  each  experiment.  The  monthly  mean  data  are 
used  for  comparison. 

4  5  Mode  Splitting 

For  computational  efficiency,  the  mode  splitting  technique  (Blumberg  and  Mellor,  1987 
is  applied  with  a  barotropic  time  step  of  25  seconds,  imsed  on  the  Courant-Friedcfichs- 
Levy  (1928;  computational  stability  (CEL)  condition  and  the  external  wave  speed;  and 
a  haroclinic  time  step  of  900  seconds,  based  on  the  CFL  condition  and  the  internal  wave 
speed. 
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4  6  Volume  Transport  Streamfunction 

To  investigate  'I'0  JES  circulation  as  a  whole,  we  integrate  the  velocity  vertically  from 
tl,!.  surface  (:  -  (»  to  the  bottom  (z  =  h) 

(I/.V)-  J  (u,v)dz  (6) 


|)u«  to  the  continuity  tin*  volume  transport  stream  function  (♦)  is  defined  by 


U  « 


<M> 

V 


V 


<MP 

dr 


and  satisfies  the  Poisson  equation 


av_  W 

dr  d\) 


(7) 

(8) 


p4)l  each  tunc  instance,  we  solve  the  two-dimensional  Pois^m  equation  ( 10)  with  the  given 
boundary  conditions,  that  is.  zero  transport  nf  the  rigid  boundaries  and  known  tran-|*irt 
at  the  open  boundaries  (see  Table  1). 


5  The  Simulated  Circulation  (Control  Run) 

5  1  General  Description 

The  simulated  surface  velocity  field  (Fig.  3)  coincides  with  the  earlier  description  of 
r lie  JF.S  circulation  presented  in  Section  2.  The  TWC  separates  at  the  I'sitshiina/Koren 
Strmt  into  western  and  eastern  branches.  Flow  through  the  western  channel  (i.e.,  EKVVC) 
closely  follows  the  Korean  coast  until  it  separates  near  38  °N  into  two  branches.  The 
eastern  branch  follows  the  SPF  to  the  western  coast  of  Snpporo  Island,  and  the  western 
branch,  moves  northward  and  forms  a  cyclonic  eddy  at  the  Eastern  Korean  Bay  (EKB). 
Flow  through  the  eastern  channel  (i.e..  .INB)  is  a  little  weaker  than  through  the  western 
channel.  File  simulated  LCC  carries  fresh  and  cold  water  along  the  llussian  coast  ami 
Utomes  the  \KCC  at  the  North  Korean  coast .  The  simulated  NKCC  meets  the  EKVVC’ 
at  about  3S>'  N  After  separation  from  the  mast,  the  NKCC  and  the  EKVVC  converge  to 
a  strong  front  that  runs  in  the  west-east  direction  across  the  basin. 

5  2  JB  Cyclonic  Gyre 

A  large-scale  cyclonic  reciiciilatinu  gyre  owr  ihe  JB  is  simulated  with  a  strong  seasonal 
variation.  This  gyre  is  easily  identified  by  tlm  volume  transport  streamfunction  (Fig.  I). 
The  JB  gyre  i»  the  strongest  and  recirculntes  8  Sv  in  the  winter.  It  weakens  and  retreats 
northward  in  tin*  spring  mid  summer.  In  the  fall,  the  cyclonic  JB  gyre  disappears  and 
weak  anticvrlomc  eddies  appears. 

5.3  LCC 

The  LCC  is  a  sonthwestwnrd  current  following  along  the  Hussian  coast  It  bifurcates  into 
twi  brandies  near  the  east  Russian  bight  (13TE,  12° N ):  the  western  branch  flows  along 
i lie  Russian -Korean  const  and  becomes  the  NKCC.  The  eastern  branch  flows  southeast 
ward,  then  turns  eastwatd  at  ll.5cN,  and  becomes  the  south  flank  of  the  JB  gyre  The 
I.CC  has  a  strong  seasonal  variation  with  the  maximum  speed  in  the  winter  (Fig.  3n) 
nnd  the  uiiniliiiini  speed  in  the  summer  (E*ig.  3c). 
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5  4  EKWC 

EKWC,  western  branch  of  TC  north  of  the  Tsushima/Korea  Strait,  acts  ns  a  western 
boundary  current  and  ha-  a  strong  seasonal  variability.  At  37'  N,  the  EKWC  v-velocity 
component,  which  is  practically  the  direction  of  the  current  on  that  position,  varies  from 
0.30  rn/s  in  winter  to  0.42  tn/s  in  summer  (Fig.  5).  The  width  of  EKWC  is  around  GO 
km  nil  year  round  However,  the  depth  is  around  1,400  tn  in  the  summer  and  800  tn  in 
the  winter. 

5.5  EKWC  -  NKCC  Confluence 

The  (northward)  overshoot  mg  of  the  EKWC  near  the  Korean  Bight  at  37  5  N  is  stronger 
in  the  winter  thnn  in  the  summer  (Fig.  3).  The  overshot  EKWC  leaves  the  Korean  coast 
and  moves  northward.  It  converges  with  the  southward  flowing,  NKCC,  at  39°N,  131/E, 
and  forms  a  current  mean  doting  toward  the  east  along  the  SPF.  We  may  call  it  the  SPF 
Current  (PFC). 

5.6  JNB 

JNB.  lie*  eastern  branch  of  TC  north  of  the  Tsushima/ Korea  Strait  at  135°E,  hns  a 
maximum  eastward  component  of  0.18  m/s  in  the  winter  (Fig.  G).  A  counter  current 
beneath  it  is  simulated  with  a  westward  component  of  0.0!)  in/.—  In  the  siunmci  the 
counter  current  remains  almost  unchanged,  but  tin*  JNB  increases  to  0.21  tn/s  (Fig  G). 

5.7  UTB  Eddy 

The  western  branch  of  the  warm  TC  (i.e.,  EKWC)  moves  northward  along  the  southern 
pait  of  Korean  coast  and  separates  from  the  coast  after  approaching  tin’  East  Korean 
Bay  (EKB).  It  keeps  its  northward  motion  until  meeting  the  NKCC  near  4()°N,  meanders 
southeastward,  and  forms  a  warm-core  anticyclonic  eddy  (Fig.  3).  The  simulated  ant  icy 
clonic  UTB  eddy  is  strongest  in  the  summer  (Fig  3c).  The  center  is  located  at  38.5°N 
131/  E.  The  size  of  the  eddy  is  around  150  km.  The  tangential  velocity  is  around  0.4  tn/s. 

6  Driving  Mechanisms 

We  analyzed  the  results  of  the  four  exj>eriinent>  to  identify  the  driving  mechanisms  for 
the  JES  circulation. 

6  1  Effects  of  Nonlinearity  (Run  1  -  Run  2) 

In  the  first  sensitivity  study,  the  nonlinear  ad  wet  ion  term-  were  removed  from  the  dy- 
namii  equations.  Otherwise,  the  same  parameters  ns  rhe  control  run  were  used.  The 
figures  and  analysis  were  performed  after  computing  the  difference  between  control  and 
non  non-linear  terms  runs.  The  differential  surface  vector  velocities  (Fig.  7)  do  not  show 
evkhuce  for  general  circulation  patterns  except  a  couple  of  mesoscale  eddies,  and  •■spe¬ 
cially  the  UTB  eddy,  which  indicates  that  rhe  non-linearity  does  not  change  the  general 
circulation  pattern  except  fot  the  UTB  anticyclonic  eddy.  However,  our  computation 
shows  that  the  non- Imparity  causes  a  noticeable  change  in  the  volume  transport. 

6.2  Wind-Induced  Circulation  (Run  1  -Run  3) 

Cotnjmml  to  the  surface  velocity  vector  map  for  the  control  rim  (Fig.  3),  the  surface 
ditferenrial  current  vector-  AV  (Fig.  8)  clearly  show  the  wind  effects  on  the  JES  circula¬ 
tion:  (a)  itriving  the  LCC  in  the  winter,  (b)  damping  the  EKWC  especially  in  the  winter, 
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u  i  generating  the  UTB  eddy  in  all  the  seasons,  and  (d)  generating  eddies  along  the  JNB. 
The  wind  has  almost  no  effect  on  the  occurrence  of  the  JNR  for  all  seasons. 

The  winter  monsoon  winds  blow  from  the  northwest  to  the  southeast  over  the  JES 
surface  and  drives  the  Ekinan  flow  in  the  upper  ocean  to  the  right  of  the  wind  direction. 
Such  a  surface  current  moves  southward  at  the  western  coast  of  the  JES  and  strengthens 
the  LCC  and  weakens  the  EKWC.  The  summer  monsoon  winds  blow  from  the  south  and 
southeast  to  the  northwest  with  much  smaller  speeds.  Thus,  the  summer  wind  effects  on 
reducing  LCC  and  strengthening  EKWC  are  quite  weak. 

The  winds  drive  mesoscale  eddies  especially  near  UTB  and  along  the  west  coast  of 
Honshu.  The  surface  UTB  eddy,  occurring  on  the  differential  velocity  vector  map  (Fig. 
SI  has  the  similar  swirl  velocity  (anticyclonic)  with  same  order  of  magnitudes  {maximum 
value  of  0.7  m/s)  as  the  control  run  (Fig.  3).  This  means  that  the  wind  effect  is  in  fact 
(jev  factor  for  generating  the  UTB  anticyclonic  eddy.  Furthermore,  there  is  no  evident 
current  along  the  west  coast  of  Honshu  on  the  AV  map  (Fig.  8},  indicating  that  wind 
forcing  is  not  a  major  factor  for  maintaining  JNB. 

6  3  Inflow/Outflow  Induced  Circulation  (Run  1  -  Run  4) 

The  third  sensitivity  study  used  the  control  run  equations  and  forcing  but  closed  all 
open  lateral  boundaries,  preventing  the  transport  of  mass,  heat  or  salinity  through  the 
Tsushima.  Tsugant,  Soya,  and  Tatar  Straits.  With  inflow  or  outflow,  the  JES  vertically 
integrated  circulation  pattern  is  more  pronounced  in  the  summer  than  in  the  winter. 
This  is  because  (control  -  no  transport)  pattern  (Fig.  9)  is  quite  similar  to  'Ifr  (control 
tuu)  (Fig-  I)  i*1  summer  than  in  the  winter.  Increased  circulation  in  the  summer  by 
the  lateral  transport  generally  leads  to  greater  horizontal  and  vertical  variabilities  of  the 
current  structure.  In  the  winter,  'I1  has  positive  values  (2  Sv),  and  fl'c  lias  negative  values 
in  JR.  I’his  suggests  that  the  lateral  transport  reduces  the  JB  cyclonic  gyre  by  2  Sv  in 
the  winter. 

7  Conclusions 

( 1 )  The  .1 ES  circulation  was  simulated  in  this  study  by  the  POM  model  under  t  he  seasonal 
surface  flux  forcing.  The  northern  JES  (i.e,.  the  Japan  Basin)  is  occupied  by  a  cyclonic 
gyre  (called  the  JB  gyre)  and  the  southern  JES  is  characterized  as  a  multi-eddy  structure. 
The  volume  ttansport  streamfunctioii  ha*  a  double-gyre  structure  with  negative  values 
(cyclonic)  in  the  northern  JES  and  positive  values  (anticyclonic)  in  the  southern  JES.  The 
simulated  .IB  gyre  is  rhe  strongest  and  recirculates  S  Sv  in  the  winter.  It  weakens  and 
retreats  northward  in  the  spring  and  summer.  In  the  fall,  the  cyclonic  JB  gyre  disappears 
and  weak  multi-anficyclonic  eddies  appear. 

(2)  The  current  systems  such  as  the  LCC,  EKWC,  NKCC,  and  JNB  are  simulated 
reasonably  well.  The  simulated  LCC  lias  a  maximum  southward  component  (0.21  m/s), 
occurring  near  the  surface  in  the  winter  with  a  width  of  100  km  and  extending  to  a  depth 
of  H|)0  in.  It  weakens  to  a  minimum  of  0.15  m/s  in  the  summer  and  fall,  and  shrinks  in 
size  to  a  width  of  60  km  and  depth  of  IIJO  tn. 

The  simulated  EKWC  varies  from  0.42  m/s  (summer)  to  0.30  rn/s  (winter).  The 
width  of  rhe  EKWC  is  around  60  km  all  year  round.  However,  the  depth  is  around  1,400 
m  in  the  summer  and  800  rn  in  the  winter.  The  (northward)  overshot  EKWC  leaves  the 
Korean  coast  moving  northward  and  converges  with  the  southward  flowing  NKCC,  at 
40* N,  130°E,  and  forms  a  current  meandering  toward  east  along  the  SPF 

The  simulated  JNB  has  a  maximum  eastward  component  of  0.24  tn/s  in  the  winter 
and  weakens  to  0.09  m/s  in  the  summer.  This  is  consistent  with  the  earlier  observational 
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si  udy.  Flu*  effect  of  coastal  geometry.  such  a*  I l>o  Nolo  Peninsula  on  the  JNB,  is  nlv. 
simulated. 

The  simulated  1TB  warm-core  ant  ky  clonic  eddy  has  a  western  branch  of  the  TWC 
•  ic.  KKWC  which  mows  northward  along  the  southern  part  of  the  Korean  coast  and 
separate*,  front  t lie  coast  after  approaching  the  East  Korean  Bay.  It  keeps  its  northward 
motion  until  meeting  the  XKCC  near  40  N.  meanders  southeastward.  and  forms  a  warm- 
core  ant tcyclonic  eddy,  This  eddy  is  strongest  in  the  summer  and  its  center  is  located  at 
38  5  N  1  -'tO '  E  The  size  of  the  eddy  is  around  150  km.  ‘Five  tangential  velocity  is  around 
0  I  m/s 

(3)  The  nonlinear  advcction  doits  not  affect  the  general  circulation  pattern  evidently 
hut  doo.  affect  the  formation  of  tin-  mesoscale  et kites,  and  es|»eciallv  the  t'TB  eddy  (all 
seasons)  and  the  JB  cyclonic  gyre  (spring). 

(4)  The  miKlel  wind  effivts  on  the  JE8  circulation  are  more  pervasive  than  those  of 
non-linear  dynamic  effects.  The  winter  winds  cause  a  strong  basin-wide  JES  cyclonic 
gyre  with  8  Sv  recirculation  in  the  northern  JES  and  the  summer  winds  drive  a  weak 
nearly  basin-wide  JES  ant  a  y  clonic  gyre  Thus,  the  wind  forcing  is  the  most  important 
fact  ;8I1%)  for  the  generation  of  the  JB  cyclonic  gyre  Besides,  the  winds  also  influence 
the  surface  circulation  such  as  driving  tlte  LCC  (winter),  damping  the  EKWC  (winter 
generating  the  ITB  eddy  (all  seasons),  and  generating  eddies  along  the  JNB  (all  seasons), 
Flic  wind  has  almost  no  effect  on  the  occurrence  of  the  JNB  for  all  seasons 

(5)  The  model  boundary-forcing  enhances  (weakens)  tlte  JES  volume  transport  in  the 
summer  (winter,.  It  has  very  weak  effects  on  the  occurrence  of  the  LCC  except  in  the 
winter,  when  the  boundary-forcing  accounts  for  30%  of  the  LCC  at  16  N.  It  weaken* 
the  JB  cyclonii  gyre  by  2  Sv  {25%)  in  the  winter  Hi  side*,  the  boundary  forcing  also 
influences  the  surface  circulation  such  as  driving  the  t1  FB  eddy  in  all  the  seasons  50% 
in  the  winter),  generating  the  EKWC  (50-100%)  in  the  winter  and  21%  in  the  summer), 
and  generating  the  JNB  and  iddies  along  the  JNB 

(6)  Future  studio*  should  concentrate  on  less  simplistic  scenarios  Realistic  lateral 
transport  should  be  included  and  the  use  of  extrapolated  climatological  winds  need*  to 
be  upgraded  to  incorporate  synoptic  winds  to  improve  realism  Finally,  the  assumption 
of  quasi-linearity  that  allowed  us  to  use  simple  differences  to  quant ify  the  effect  of  ex¬ 
ternal  forcing  needs  to  be  rigorously  tested,  It  is  very  important  to  develop  a  thorough 
methodology  to  perform  sensitivity  studies  under  the  highly  non-linear  conditions  that 
may  exist  in  the  littoral  environment, 
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E'lg.  1.  Geography  and  isobaths  showing  the  bathymetry  (m)  of  the  Japan/East 
Sea  (JES). 
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Fig.  2.  Schematic  map  of  surface  current  systems  (after  Tomczak  and 
Godfrey,  1994). 
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Fig.  3.  Simulated  surface  circulation  from  the  control  run  for  (a)  January, 
(b)  April,  (c)  July  and  (d)  October. 
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Fig.  4.  Simulated  volume  transport  (Sv)  from  the  control  run  for  (a) 
January,  (b)  April,  (c)  July  and  (d)  October. 
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Fig.  7.  Surface  circulation  caused  by  the 
(b)  April,  (c)  July  and  (d)  October. 
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Fig.  8.  Surface  circulation  driven  by  the  winds:  (a)  January,  (b)  April,  (c) 
July  and  (d)  October. 
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Fig.  9.  Volume  transport  driven  by  the  lateral  forcing:  (a)  January,  (b) 
April,  (c)  July  and  (d)  October. 
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